Antiferromagnetic Critical Fluctuations in BaFe 2 As 2 
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Magnetic correlations near the magneto-structural phase transition in the bilayer iron pnictide 
parent compound, BaFe2As2, are measured. In close proximity to the antiferromagnetic phase 
transition in BaFe2As2, a crossover to three dimensional critical behavior is anticipated and has been 
preliminarily observed. Here we report complementary measurements of two-dimensional magnetic 
fluctuations over a broad temperature range about Tjv. The potential role of two-dimensional critical 
fluctuations in the magnetic phase behavior of BaFe2As2 and their evolution near the anticipated 
crossover to three dimensional critical behavior and long-range order are discussed. 



PACS numbers: 74.70.Dd, 75.25.Z, 75.50.Ee, 75.40.Cx 
I. INTRODUCTION 

Reduced dimensionality constitutes a key element of 
the magnetism inherent in the lamellar copper oxides 
and is also thought to play a key role within the mech- 
anism responsible for the formation of the high tem- 
perature superconducting (high-T c ) phase 1 . The recent 
discovery of superconductivity at temperatures as high 
as T c = 55K within an alternate class of iron based 
superconductors 2-5 has stimulated a host of compar- 
isons between the ferrous parent compounds of these 
new superconductors and the Mott-insulating parent 
systems of the cuprates. In particular, antiferromag- 
netic order exhibiting energetic and anisotropic magnetic 
exchange interactions is observed within both classes 
of materials 1, 6-8 . While the spin ordering in high-T c 
cuprates such as La2Cu04 has been demonstrated to 
arise from strong two-dimensional (2D) spin fluctuations 
consistent with a 2D-Hcisenberg model with weak inter- 
layer coupling 6 ' 10 , the buildup of correlations yielding the 
antiferromagnetic parent phase in the iron based super- 
conductors remains relatively unexplored. 

Specifically, within the bilayer iron pnictide variant 
BaFe2As2 (Ba-122) studies of both pure and doped sam- 
ples have suggested a second order or weakly first order 
phase transition at Tjv = 136K n ~ 15 . Critical fluctua- 
tions in proximity to the onset of long-range magnetic 
order have long been known to yield important insight 
into the fundamental spin and spatial symmetries driv- 
ing the phase formation. It is therefore important to ex- 
plore the nature of the fluctuations driving the magnetic 
phase transition in Ba-122 as a means of understanding 
the spin behavior in this iron pnictide parent system as it 



passes through its simultaneous magnetic and structural 
phase transitions. 

BaFc 2 As 2 develops long range antiferromagnetic or- 
der at T= 136K concomitant to a structural distortion 
from tetragonal to orthorhombic symmetry 16 . The bi- 
quadratic free energy coupling term between the struc- 
tural and magnetic order parameters (implied by their 
identical phase behaviors) suggests an unorthodox cou- 
pling scheme rather than a simple linear coupling of 
the magnetic order parameter squared to the lattice 
through strain 11,17 ' 18 . Models such as fluctuating mag- 
netic domains locally breaking time reversal symmetry 
prior to the transition to long-range order 19 or the sup- 
position of a weakly first order nature to the structural 
phase transition may explain this peculiar coupling term; 
however other symmetries such as higher energy orbital 
ordering 20 ' 21 or valley density waves 22 must be consid- 
ered as well. 

In this article, we present a neutron scattering study 
of 2D magnetic fluctuations above and below the antifer- 
romagnetic phase transition in BaFe2As2- Given the 2D- 
Ising-like exponent of the magnetic phase transition 11,23 , 
2D-Heisenberg fluctuations are expected at temperatures 
far above Tjy. As the system is cooled sufficiently close to 
Tjv, the spin system must crossover to three-dimensional 
behavior, and the longitudinal component of these two- 
dimensional fluctuations is expected to condense into the 
three-dimensional Bragg positions for the antiferromag- 
netic ordered phase. Transverse spin fluctuations how- 
ever are expected to remain noncritical and simply evolve 
continuously into the spin wave component of the spec- 
trum as the system is tuned across Tjv- Here we present 
new measurements of the energy-integrated, two dimen- 
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FIG. 1. (a) Schematic showing the scattering configuration 
employed to measure the energy integrated, instantaneous 
correlation length with c || kf (b) Two-dimensional spin struc- 
ture for BaFe2As2 showing spin orientations within the Fe- 
plane only. Critical fluctuations measured within the spec- 
trometer configuration in panel (a) at non-integer L necessar- 
ily arise from uncorrelated planes with this in-plane ordering 
pattern. 



sional fluctuations in Ba-122 that are consistent with 
this picture, and these results complement our previ- 
ous observation of 3D critical fluctuations persisting until 
T« 150K (1 - ^- =t« 0.10) in the same sample 11 . Our 
study supports trie notion that the magnetic phase tran- 
sition in the BaFe2As2 system is driven fundamentally 
by magnetic fluctuations of two dimensional character. 



II. EXPERIMENTAL 

A 600mg crystal of BaFe2As 2 was grown out of an 
FeAs flux with a ratio Ba:Fe:As of 1:4:4 and possessed 
a mosaicity of less than 0.3° FWHM. The crystal was 
mounted in the [H, 0, L] plane (orthorhombic notation 
a = b = 5.6202A, c = 13.0323A at T= 300K) in- 
side of a closed cycle refrigerator. Positions in recipro- 
cal space are given in reduced lattice notation where Q 
[r.l.u.] =Q[A _1 ]-(^pft + J^k + ^A). Experiments were 
performed on the BT-7 triple-axis spectrometer at the 
NIST Center for Neutron Research at the National In- 
stitute of Standards and Technology and on the HB-1 
triple-axis spectrometer at the High Flux Isotope Reac- 
tor at Oak Ridge National Lab. Experiments on BT- 
7 were performed with an incident neutron energy of 
E< = 14.7meV and collimations of open-50'-sample-50'- 
detector. Experiments on HB-1 were performed with an 
incident energy of E< = 13.5meV and collimations of 40'- 
40'-sample-60'-240'-detector. In all experiments, two PG 
filters were placed before the sample in order to elimi- 
nate higher order contamination. The (0, 0, 2) reflec- 
tions of pyrolitic graphite (PG) crystals were used as the 
incident energy on vertically focusing monochromaters. 
Uncertainties presented in the data are statistical in ori- 
gin and represent one standard deviation. Experiments 
were performed with the analyzer crystal removed so that 
two-axis energy integration along the k/ side could be 



performed. 



III. MAGNETIC SCATTERING IN BaFe 2 As 2 

In order to isolate the inherently 2D component of the 
spin scattering in Ba-122 resulting from magnetism corre- 
lated only within the Fe-planes, the experimental geom- 
etry was chosen such that the c-axis was parallel to the 
final scattering wave vector, k/ 24,25 . The analyzer was 
removed thus ensuring that any two-dimensional diffuse 
scattering along L at finite energy transfers maintains the 
specified in-plane H- value (Fig. 1). When integrating 
over various energy transfers at a fixed two-theta value, 
the corresponding momentum transfers probed by the 
neutron must also vary; however through orienting c || k/ 
the component of momentum transfer that varies with fi- 
nal energy is constrained to be along only the L-direction 
in momentum space. For scattering that is diffuse along 
L, this technique is uniquely powerful in providing an ef- 
ficient tool for the measurement of an energy integrated 
spin response. This allows the instantaneous correlation 
function for any 2D component of magnetic scattering 
within the Fe-planes to be measured through scans along 
H with fixed k^-dependent L-values. Furthermore, our 
geometry is also an efficient method for probing the pres- 
ence of diffuse, quasi-2D scattering through searching 
for magnetic signal at nominally disallowed 3D magnetic 
Bragg positions. For instance, when probing the in-plane 
magnetic zone center at H= 1 with Ej = 14.7 meV, this 
translates to a geometry where the c-axis is parallel to 
kf at an L-value of L= 0.51. 

Energy transfers probed by neutron energy loss are 
integrated up Ej while all energy transfers probed by 
neutron energy gain are integrated. The most relevent 
range of integrated L values, assuming neutron energy 
transfers in the window from —11.7 < AE < 14.7 meV, 
can be estimated as —5.01 <L>i nt < 2.39 (r.l.u.). This 
range accounts for the most significant thermally acti- 
vated fluctuations at T^r. Thus, the presence of cor- 
related magnetic scattering at the (1, 0, 0.51) position 
in this geometry provides qualitative evidence of 2D dif- 
fuse scattering integrated within the detector rather than 
dispersive 3D excitations that would cut through only a 
small fraction of our experimental window and be lost in 
the larger background incurred via this technique. 

In Fig. 2 (a), we have plotted the raw data for scans 
along the [H, 0, 0] direction with L=0.51 (this satisfies 
c|| kf for Ej = 14.7meV). Data taken at T=Tjv are over- 
plotted with data collected far below the transition tem- 
perature at T=80K. From these data, it is immediately 
evident that at the transition a substantial amount of 
scattering appears centered at the Qi n - P i an e — (1> 0) po- 
sition. Subsequent H-scans for T< 100K show no further 
changes in intensity at Qin-piane — (1,0) upon cooling, 
indicating that the 2D fluctuations have disappeared by 
this temperature. Both H- and L-scans below these tem- 
peratures (blue symbols) therefore serve as a useful back- 
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ground for the removal of nonmagnetic scattering and for 
the isolation of critical scattering near 7jv. Similarly, the 
small peak at H= 1.25 in Fig. 2 (a) is a temperature 
independent nuclear peak from another small crystallite 
with the sample. The background subtracted data at 
T= 136K in Fig. 2 (c) show a clear peak centered at the 
Qin-piane — (1,0) position indicative of scattering aris- 
ing from quasi-elastic 2D, in-plane fluctuations at the 
AF transition. In order to prove that the signal observed 
at Qin-piane = (1,0) is indeed two-dimensional, scans 
were performed along L, and plotted in Figs. 2 (b) and 
(d) showing corresponding raw and subtracted L-scans 
respectively. The subtracted scan in Fig. 2 (d) shows 
the critical scattering peaks at the two-axis integration 
condition of c || kf at L= 0.51 — as expected for 2D dif- 
fuse scattering at the Qin-piane = (1,0) position. This 
broad peak, extending over half of the Brillouin zone in 
L, sharply contrasts with the relatively sharp peak ob- 
served in H-scans and is indicative of 2D fluctuations 
arising from spins associated with the antiferromagnetic 
phase. However, it is important to note that the width of 
the peak in L is predominantly controlled by the relative 
orientation of c with respect to kf (ie. the degree of inte- 
gration performed rather than the correlations between 
the planes). The intrinsic spectrometer resolution is a 
neglible component of the L-scan's peak width as shown 
in Fig. 2. 

2D short-range spin correlations persist far above Tjv 
at the AF ordering wave vector and eventually vanish 
by sa 250K. Figure 3 (a) shows background-subtracted 
H-scans through the (1,0) position at 136K and 300K, 
revealing the absence of the 2D scattering signal at high 
temperatures. A gradual increase in the overall back- 
ground intensity was observed above 150K and is re- 
flected by the large uniform offset in the 300K subtracted 
data. 

In investigating the detailed temperature dependence 
of the critical scattering signal, the peak intensity of 
the scattering at (1,0) was measured as the sample was 
slowly warmed through Tjv. The resulting intensity as a 
function of temperature is plotted in Fig. 3 (b) where the 
intensity of the 2D scattering sharply increases upon ap- 
proaching Tjv from below and gradually decreases as the 
system is warmed away from the transition temperature. 
In order to compensate for the increase in background 
above 150K, the nonmagnetic background contribution 
to the scattering was measured by rotating the crystal 
15° away from the (1, 0)i n ~ p iane position, and it was 
subsequently removed from the data shown in Fig. 3 
(b). By 250K the peak has broadened substantially and 
is beyond the experiment's ability to resolve it above the 
background. 

The intensity of the 2D spin scattering peaks in close 
proximity to 136K near the onset of 3D AF order and is 
consistent with the scattering intensities at finite energy 
transfer reported in CaFe2As2 (Ca-122) 26 . To provide 
a better reference for this, the purely elastic (AE= 0) 
intensity of the Q= (1,0,3) position is overplotted with 
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FIG. 2. (a) H-scans through the energy integration position 
(1, 0, 0.51) both below (blue square symbols) and at T/v 
(black circle symbols), (b) Raw data from L-scans through 
(1, 0, 0.51) below and at Tjv respectively. Background- 
subtracted data for both H-scans and L-scans through (1, 0, 
0.51) are also plotted in panels (c) and (d) respectively. Verti- 
cal dashed lines in panels (c) and (d) illustrate the full width 
at half maximum (FWHM) of the spectrometer resolution 
function projected along the H and L directions respectively. 



the 2D scattering intensity in Fig. 3 (b). The sharp de- 
crease in 2D scattering intensity below Tjv likely reflects 
the spin gap opening below Tjy and the subsequently re- 
duced scattering intensity as it continuously approaches 
its full A = lOmeV value 13 . The buildup in scatter- 
ing intensity in approaching Tjv from above reflects the 
buildup of in-plane correlations, both longitudinal and 
transverse components, as the system is cooled toward 
the phase transition. The dominant contribution to this 
signal should originate from the transverse, spin wave, 
component of the 2D scattering, and this will be dis- 
cussed in more detail in the Discussion section of this 
paper. 

We next investigated the manner through which the 
correlation length of the 2D fluctuations evolved as the 
sample was cooled toward Tjv- In principle, this rep- 
resents an approximation of the energy-integrated, in- 
stantaneous correlation function; however the energetic 
spin excitations in BaFe2As2 preclude a total integra- 
tion. Instead, only the most relevant, low-energy, fluc- 
tuations are probed with the incident energies used. 
The background-subtracted data were fit to a two- 
dimensional Lorentzian profile convolved with the instru- 
mental resolution. Within this convolution, the elastic 
resolution of the spectrometer was convolved into the fit 
peaks, which effectively underestimates the total resolu- 
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FIG. 3. Critical scattering signal in BaFe2As2. (a) H-scans 
showing background subtracted data at 136K (black circles) 
and 300K (red squares) collected on BT-7. (b) BT-7 data 
showing the background subtracted energy integrated inten- 
sities of 2D fluctuations plotted as a function of temperature 
(red circles). The uncertainty for the green circle shows the 
uncertainty in a vanishing signal at 300K. The vertical dashed 
line denotes Tjv, and the solid black line through the data is 
a guide to the eye. Elastic data collected at Q=(l, 0, 3) in 
a triple-axis configuration on HB-1 are also plotted as black 
triangles and show the development of the magnetic order 
parameter in the sample. 



tion contribution to the scattering profile. The differing 
resolution functions for dynamic correlations from finite 
energy transfers must be accounted for; however without 
knowing the precise spectral distribution of the local spin 
susceptibility for each corresponding temperature, this is 
not possible. Instead, the correlation lengths plotted rep- 
resent the minimum spin-spin correlation length at each 
temperature. The results of fits at various temperatures 
are shown in Fig. 4. At high temperatures far from Tjy 
the correlations are short-ranged, and as the system is 
cooled toward T^r the correlation lengths begin to in- 
crease. Correlation lengths smoothly increase upon de- 
creasing temperature and seemingly saturate around Tjv. 
The absence of a true divergence to long-range order at 
the phase transition is expected given that the critical, 
longitudinal, component of the scattering necessary con- 
denses into the 3D Bragg position as the system crosses 
over from 2D to 3D regimes above Tpf. Hence, near Tjv 




Temperature (K) 

FIG. 4. 2D, in-plane, spin correlation lengths measured in the 
two-axis configuration both above and below Tjv- The verti- 
cal dashed line denotes Tjv. Both E; = 13.5 and 14.7 meV 
were used and scans conducted along H centered at the corre- 
sponding (1, 0, L) integration positions. Data were collected 
on both the HB-1 and BT-7 spectrometers, and the solid line 
is a guide to the eye. 



we are measuring the noncritical, transverse correlation 
length. 



IV. DISCUSSION 

The two-dimensional critical scattering observed in 
BaFe2As2 far above T^r is consistent with the 2D Heisen- 
berg spin excitations anticipated at temperatures far 
above the 2D Ising-like magnetic order parameter in this 
system. At these high temperatures, the scattering rep- 
resents a mixture of both longitudinal and transverse 
2D critical fluctuations. Projecting the components of 
the measured susceptibility for an untwinned crystal into 
vectors perpendicular and parallel to the ordered mo- 
ment direction (?) = (1,0,0)), the measured spin sus- 
ceptibility can be expressed as the linear combination 
= - sin 2 {9)) + xHsin 2 (6)) where 9 is the an- 

gle between the Q-vector and the a-axis and x 1 ' an d 
are projections of the spin susceptibility within the 
(b, c)-plane and along the a-axis respectively 25 . For a 
twinned crystal at the Q= (1,0,0.51) position this be- 
comes |^ = lA77x ± +0.523x" where approximately 5% 
of the measured signal is comprised of longitudinal fluc- 
tuations, x ■ 

Upon cooling, the longitudinal component of the crit- 
ical scattering is expected to condense into the 3D AF 
Bragg positions at the dimensional crossover T30 rj 
150K inferred from our previous studies' . The remain- 
ing transverse component represents the spin wave contri- 
bution to x" which is nondivergent at Tn- This picture 
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is consistent with the observed behavior of the 2D fluc- 
tuations in both Figs. 3 and 4. Due to the minority com- 
ponent of x" in the measured 2D signal and the known 
crossover in the dimensionality of this system away from 
T^, the correlation lengths plotted in Fig. 4 reflect the 
evolution of the spin-spin correlations of in-plane spin 
wave excitations. Below Tjv within the resolution of our 
experiments, we could detect no further change in the 2D 
correlation lengths; however a rapid decrease in scatter- 
ing signal prevented measurements far below TV- 

The buildup in low energy spectral weight upon ap- 
proaching Tjv from above is consistent with the spin be- 
havior observed within isostructural Ca-122 26 and with 
the expected divergence of the staggered spin susceptibil- 
ity upon approaching T N . Upon cooling from high tem- 
peratures, 2D critical fluctuations build up within the 
longitudinal component which are subsequently inter- 
rupted at T 3 £>, and our data suggest that the remain- 
ing transverse fluctuations continue to build smoothly 
toward TV. This conjecture is based on our expecta- 
tion that the nondivergent behavior at T^ is suggestive 
of remnant transverse fluctuations. Below Tjv, the inten- 
sity is therefore representative of only the transverse crit- 
ical fluctuations that are subsequently diminished upon 
further cooling. The suppression of quasi-2D scattering 
upon cooling below T N is primarily due to the opening of 
a spin anisotropy gap within the ordered phase 13 whose 
value increases rapidly away from the 3D zone center with 
the buildup of intcrplane correlations below Tn and the 
resulting dispersion of spin excitations along L. 

Another issue worth commenting on is the lack of any 
detectable anomaly in the quasi-2D signal in Fig. 3 at 
T= 150K upon approaching T^ from above. Naively, 
one may expect to resolve some shift in the intensity of 
the 2D fluctuations as the longitudinal component of the 
signal condenses into the 3D AF ordering wave vectors 
near T 3 £>. Dimensional crossovers however typically oc- 
cur over the span of a decade in reduced temperature 
rendering changes resolvable primarily through shifts in 
the critical exponents of the phase transition. Therefore, 
our current data preclude any definitive statement on 
the presence or absence of a crossover primarily due to 
insufficient statistics, although future experiments with 
increased sample volume will help resolve this issue. 

In attempting to analyze the data below T^ , one must 
consider the rapidly evolving dispersion as the system is 



cooled through both the magnetic and structural phase 
transitions at T= 136K. At present, there has been no 
detailed study of the evolution of magnetic exchange cou- 
plings where the nearest neighbor Jia~Jib above T^ 
transitions into a highly anisotropic antifcrromagnetic 
Jia and ferromagnetic Jib below Tjy 7 . Until the na- 
ture of this transition is known in detail, the complexity 
of the magnetic phase transition precludes a more in- 
depth analysis of the critical behavior below Tjv- The 
small sample volume of the Ba-122 crystal in the cur- 
rent study is the limiting factor rather than an intrinsi- 
cally weak critical scattering signal 24 , and future studies 
with a greatly increased sample volume will allow a more 
comprehensive evaluation of the critical behavior in this 
material. 



V. CONCLUSIONS 

We have investigated the behavior of 2D spin fluctua- 
tions about the magnetic phase transition in the Ba-122 
system. From this initial study, we have resolved a mix- 
ture of both transverse and longitudinal critical fluctua- 
tions whose onset is coupled to appearance of 3D AF at 
Tjv- This suggests that 2D fluctuations play an impor- 
tant role in driving the formation of the magnetic phase 
within this system; however the interrupted divergence 
in the instantaneous correlation lengths and crossover 
to 3D spin behavior within this 2D magnetic scattering 
require further investigation. Future experiments with 
larger sample volume and tighter experimental collima- 
tion will allow a more quantitative assessment of energy 
integrated, spin-spin correlation lengths close to T N . 
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